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Abstract 

This paper is a short review on the application of continuos-time random 
walks to Econophysics in the last five years. 

1 Introduction 

Recently, there has been an increasing interest on the statistical properties 
of high-frequency financial data related to market microstructural properties 
El EH El El El • High-frequency econometrics is now well established after 
research on autoregressive conditional duration models El El QUI ■ 

In high-frequency financial data not only returns but also waiting times 
between consecutive trades are random variables |llj . This remark is present 
in a paper by Lo and McKinlay ^2Jj but it can be traced at least to pa- 
pers on the application of compound Poisson processes and subordinated 
stochastic processes to finance. Models of tick-by-tick financial data based 
on compound Poisson processes can also be found in the following references: 

E3E1G21. 

Compound Poisson processes are an instance of continuous-time random 
walks (CTRWs) ^H]. The application of CTRW to economical problems dates 
back, at least, to the 1980s. In 1984, Rudolf Hilfer discussed the application 
of stochastic processes to operational planning, and used CTRWs as tools for 
sale forecasts U^]. The revisited and augmented CTRW formalism has been 
applied to high-frequency price dynamics in financial markets by our research 
group since 2000, in a series of three papers |20l 1211 12*2*] . Other scholars have 
recently used this formalism |23l 1241 12*5] . However, already in 1903, the PhD 
thesis of Filip Lundberg presented a model for ruin theory of insurance compa- 
nies, which was later developed by Cramer |26H27| . The underlying stochastic 
process of the Lundberg-Cramer model is another example of compound Pois- 
son process and thus also of CTRW. 
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Among other issues, we have studied the independence between log-returns 
and waiting times for stocks traded at the New York Stock Exchange in Oc- 
tober 1999. For instance, according to a contingency-table analysis performed 
on General Electric (GE) prices, the null hypothesis of independence can be 
rejected with a significance level of 1 % [2SI - We have also discussed the 
anomalous non-exponential behaviour of the unconditional waiting-time dis- 
tribution between tick-by-tick trades both for future markets [21] and for stock 
markets 28 , 29 . Different waiting-time scales have been investigated in differ- 
ent markets by various authors. All these empirical analyses corroborate the 
waiting-time anomalous behaviour. A study on the waiting times in a contem- 
porary FOREX exchange and in the XlXth century Irish stock market was 
presented by Sabatelli et al. [S3j • They were able to fit the Irish data by means 
of a Mittag-Lefficr function as we did before in a paper on the waiting-time 
marginal distribution in the German-bund future market |21| . Kyungsik Kim 
and Seong-Min Yoon studied the tick dynamical behavior of the bond futures 
in Korean Futures Exchange (KOFEX) market and found that the survival 
probability displays a stretched-exponential form (23- Finally, Ivanov et al. 
33 confirmed that a stretched exponential fits well the survival distribution 
for NYSE stocks as we suggested in [2E|- Moreover, just to stress the relevance 
of non-exponential waiting times, a power-law distribution has been recently 
detected by T. Kaizoji and M. Kaizoji in analyzing the calm time interval 
of price changes in the Japanese market [22] • We have offered a possible ex- 
planation of the anomalous waiting-time behaviour in terms of daily variable 
activity 

The aforementioned empirical results are important as market microstruc- 
tural models should be able to reproduce such a non-exponential behaviour 
of waiting-time distributions in order to be realistic. However, the rest of this 
paper focuses on the theory and is divided as follows: in Section 2, CTRW 
theory is presented as applied to finance. Finally, in Sec. 3, a summary of 
main results is presented together with a discussion on the direction of future 
research. 



2 Theory 

Random walks have been used in finance since the seminal thesis of Bachelier 
[231 , a work completed at the end of the XlXth century, more than a hundred 
years ago. After a rather long period in which the ideas of Bachelier were 
neglected, they were further developed until recent times |35l l3l)| . 

Our approach to random walks in finance is related to that of Clark [T^j 
and to the introductory part of Parkinson's paper |37j . It is a purely phe- 
nomenological approach. There is no assumption on the rationality or the 
behaviour of trading agents and it is not necessary to assume the validity of 
the efficient market hypothesis However, as briefly discussed above, 

even in the absence of a microfoundation, a phcnomcnological model can still 
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be useful to corroborate or falsify the consequences of behavioural or other 
assumptions on markets. The model itself can be corroborated or falsified by 
empirical data. 

In order to model tick-by-tick data, we use the so-called continuous-time 
random walk (CTRW), where time intervals between successive steps are ran- 
dom variables, as discussed by Montroll and Weiss ^Hj- In physics, CTRWs 
have been introduced as models of diffusion with instantaneous jumps from 
one position to the next. For this reason they can be used as models of price 
dynamics as well. 

Let S(t) denote the price of an asset at time t. In a real market with a 
double-auction mechanism, prices are fixed when buy orders are matched with 
sell orders and a transaction (trade) occurs. It is more convenient to refer to 
returns rather than prices. For this reason, we shall take into account the 
variable x(t) = logS'(t): the logarithm of the price. For a small price variation 
AS = S(ti + i) — S(ti), the return r — AS/S(ti) and the logarithmic return 
flog = log[S(ti+i)/ S(U)] virtually coincide. 

CTRWs are essentially point processes with reward 00] • The point process 
is characterized by a sequence of independent identically distributed (i.i.d.) 
positive random variables T{, which can be interpreted as waiting times be- 
tween two consecutive events: 

n 

t n = h + '^Ti\ t n — t n -x = T n ; n= 1,2,3,...; t = 0. (1) 

i=l 

The rewards are (i.i.d.) not necessarily positive random variables: In the 
usual physical intepretation, the £jS represent the jumps of a diffusing par- 
ticle, and they can be n-dimensional vectors. Here, only the 1-dimensional 
case is studied, but the extension of many results to the n-dimensional case 
is straightforward. The position x of the particle at time t is given by the 
following random sum (with N(t) — max{n : t n < t} and a;(0) = 0): 

N(t) 

ar(t) = X> (2) 

i=l 

In the financial interpretation outlined above, the £j's have the meaning of log- 
returns, whereas the positions or rewards x(t) represent log-prices at time t. 
Indeed, the time series {x(U)} is characterised by r), the joint probability 
density of log-returns £j = x(ti + i) — x(ti) and of waiting times r, = — U. 
The joint density satisfies the normalization condition J f d^drip(^,T) = 1. It 
must be again remarked that both £j and Tj are assumed to be independent 
and identically distributed (i.i.d.) random variables. This strong assumption 
is useful to derive limit theorems for the stochastic processes described by 
CTRWs. However, in financial time series, the presence of volatility clustering, 
as well as correlations between waiting times do falsify the i.i.d hypothesis. 
The reader interested in a review on correlated random variables in finance is 
referred to chapter 8 in McCauley's recent book |41| . 
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In general, log-returns and waiting times are not independent from each 
other [2HI- By probabilistic arguments (see |18l 1211 I32| L one can derive the 
following integral equation that gives the probability density, p(x,t), for the 
particle of being in position x at time t, conditioned by the fact that it was 
in position x = at time t = 0: 

rt r + co 

p(x,t) = 5(x)W(t)+ / <p(x - x',t-t')p(x',t')dt' dx', (3) 



where 6(x) is Dirac's delta function and <?(r) is the so-called survival function. 
\P(t) is related to the marginal waiting-time probability density iP(t). The two 
marginal densities ip{r) and A(£) are: 

r+ao 



A(0 = / <p{Z,T)&r, (4) 
Jo 

and the survival function S'(r) is: 

i>T poo 

!?(r) = l-y tp{T')dr' = J 4>{T')dT'. (5) 

Both the two marginal densities and the survival function can be empirically 
derived from tick-by-tick financial data in a direct way. 

The integral equation, eq. can be solved in the Laplace- Fourier domain. 
The Laplace transform, g(s) of a (generalized) function g{t) is defined as: 



g(s)= / dte~ st g(t), (6) 
Jo 

whereas the Fourier transform of a (generalized) function f(x) is defined as: 

r+oo 

/(«)= / dxe iK *f(x). (7) 



A generalized function is a distribution (like Dirac's 5) in the sense of S. L. 
Sobolev and L. Schwartz l.'i . 
One gets: 

%K,8)=#{S) J- -, (8) 

1 - <p(K,S) 

or, in terms of the density i^(t): 

p(K,s) = f^l -, (9) 

s l-<p(K,S) 

as, from eq. (0, one has: 
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= i-^ 

s 

In order to obtain p(x,t), it is then necessary to invert its Laplace- Fourier 
transform p(n,s). As we shall see in the next subsection, for log-returns in- 
dependent from waiting times, it is possible to derive a series solution to the 
integral equation @- 



2.1 Limit Theorems: The Uncoupled Case 

In a recent paper, Gorenflo, Mainardi and the present author have discussed 
the case in which log-returns and waiting times are independent [|2 ■ It i s 
the so-called uncoupled case, when it is possible to write the joint probability 
density of log-returns and waiting times as the product of the two marginal 
densities: 

v(e,r)=A(OV(r) (11) 

with the normalization conditions f d£\(£) — 1 and / dTip(r) = 1. 
In this case the integral master equation for p(x, t) becomes: 



\(x - x')p(x',t')dx' 



dt' (12) 



This equation has a well known general explicit solution in terms of P(n,t), 
the probability of n jumps occurring up to time t, and of the n-fold convolution 
of the jump density, X n (x): 

/+oo r+oo 
... deiA(a:-&,_i)...A(6). (13) 

-oo J — oo 

Indeed, P(n, t) is given by: 

P(n,t) = f i> n (t-T)V{T)dT (14) 
Jo 

where tp n {T) is the n-fold convolution of the waiting-time density: 

-0n(l~) = / .../ dT n -l ...dTll(>(T-T n -l)...ll>(Tl). (15) 



The n-fold convolutions defined above are probability density functions for 
the sum of n independent variables. 

Using the Laplace-Fourier method and recalling the properties of Laplace- 
Fourier transforms of convolutions, one gets the following solution of the in- 
tegral equation |31 1521 1351 : 

OO 

p[x,t) = Y,P{n,t)X n {x) (16) 

n=0 
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Eq. i|16|) can also be used as the starting point to derive eq. (|12fl via the 
transforms of Fourier and Laplace, as it describes a jump process subordinated 
to a renewal process [T11IT7] . 

Let us now consider the following pseudodifferential equation, giving rise 
to anomalous relaxation and power-law tails in the waiting-time probability: 

dP 

—9(t) = -W(t), t>0, 0</3<l; 9(0+) = 1, (17) 

where the operator dP /dt 13 is the Caputo fractional derivative, related to 
the Riemann-Liouville fractional derivative. For a sufficiently well-behaved 
function fit), the Caputo derivative is defined by the following equation, for 
<{3 < 1: 

# ^ 1 d r f(r) ,_ t-P _ + . 



WW = W^fT) dt J jf^ry dT - T^W)^ ^ 

and reduces to the ordinary first derivative for (3 = 1. The Laplace transform 
of the Caputo derivative of a function fit) is: 

£ (!?/(*); s ) = s "f( s ) ~ ^V(o+). (i9) 

If eq. (|19f) is applied to the Cauchy problem of eq. l|17f) . one gets: 

*w = (20) 

Eq. (|2(J[1 can be inverted, giving the solution of eq. Q17[) in terms of the Mittag- 
Leffler function of parameter (3 |481 0H| : 

9{t) = £,3(-A (21) 

defined by the following power series in the complex plane: 

oo n 

E iz) := V — r . (22) 

n—0 v ' 

The Mittag-Leffler function is a possible model for a fat-tailed survival func- 
tion. For (3 — 1, the Mittag-Leffler function coincides with the ordinary ex- 
ponential function. For small r, the Mittag-Leffler survival function coincides 
with the stretched exponential: 

V(t) = E p (-t?) ~ 1 - ~ eM-r /r(/3 + 1)}, < r « 1, (23) 

whereas for large r, it has the asymptotic representation: 
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, s sin(/37r) r(B) , . 

#( T ) 0</3<l, t^oo. (24) 

Accordingly, for small r, the probability density function of waiting times 
V>(t) = —d^{r)/dr behaves as: 

W = - < r « 1, (25) 

and the asymptotic representation is: 

, sin(/37r) r(8 + l) , x 

V(r) W-J. ^ \ o < 3 < 1, r -f oo. (26) 

The Mittag-Lefncr function is important as, without passage to the diffu- 
sion limit, it leads to a time-fractional master equation, just by insertion into 
the CTRW integral equation. This fact was discovered and made explicit for 
the first time in 1995 by Hilfer and Anton [HOj- Therefore, this special type 
of waiting-time law (with its particular properties of being singular at zero, 
completely monotonic and long-tailed) may be best suited for approximate 
CTRW Monte Carlo simulations of fractional diffusion. 

For processes with survival function given by the Mittag-Leffler function, 
the solution of the master equation can be explicitly written: 

P(M) = £— Ef\- t P)\ n {x), (27) 



ni 

n=0 



where: 



The Fourier transform of eq. (|27l) is the characteristic function of p(x, t) 
and is given by: 

p(K,t)=E^{X(K)-l)}. (28) 
If log-returns and waiting times are scaled according to: 

x n (h) = h£x + h& + ■ ■ ■ + h£ n , (29) 

and: 

t n (r) = tt\ + rr 2 + . . . + rr n , (30) 
the scaled characteristic function becomes: 



Ph,r(K,t) = E 



73 



-j(A(ft/s) - 1; 



(31) 



Now, if we assume the following asymptotic behaviours for vanishing h and r: 
\{hn) ~ l-h a \n\ a ; < a < 2, (32) 
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and 

h a 

lim —- = 1, (33) 

h.r^O rP 

we get that: 

lim p h , r (K,t) =u(K,t) =E [-tP\K\ a }. (34) 



The Laplace transform of eq. I|34|) is: 

u(k, s) = t— p a- (35) 



Therefore, the well-scaled limit of the CTRW characteristic function coincides 
with the Green function of the following pseudodifferential fractional diffusion 
equation: 

\K\ a i{K, s) + s^(k, s) ^s?- 1 , (36) 

with u(x, i) given by: 

u{x,t) = ^W a>p ^), (37) 



where W a> p(u) is given by: 



1 



+ 00 



W a ,p(i*) = 7T / d K e- lKU E f3 (-\ K \ a ), (38) 

the inverse Fourier transform of a Mittag-Leffler function |52 E3 I5H ESI • 
For /3 = 1 and a = 2, the fractional diffusion equation reduces to the 
ordinary diffusion equation and the function W2,i(u) becomes the Gaussian 
probability density function evolving in time with a variance a 2 = 2t. In the 
general case (0 < [3 < 1 and < a < 2), the function W a ,p{u) is still a 
probability density evolving in time, and it belongs to the class of Fox ii-type 
functions that can be expressed in terms of a Mellin-Barnes integral as shown 
in details in ref. |52j . 

The scaling equation, eq. (|33|l . can be written in the following form: 

h ~ r' 9 /". (39) 

If P = 1 and a — 2, one recognizes the scaling relation typical of Brownian 
motion (or the Wiener process). 

In the passage to the limit outlined above, p r h {n, s) and u(k, s) are asymp- 
totically equivalent in the Laplace-Fourier domain. Then, the asymptotic 
equivalence in the space-time domain between the master equation and the 
fractional diffusion equation is due to the continuity theorem for sequences 
of characteristic functions, after the application of the analogous theorem 
for sequences of Laplace transforms 56 . Therefore, there is convergence in 
law or weak convergence for the corresponding probability distributions and 
densities. Here, weak convergence means that the Laplace transform and/or 
Fourier transform (characteristic function) of the probability density function 
are pointwise convergent (see ref. |56|V 
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2.2 Limit Theorems: The Coupled Case 

The diffusive limit in the coupled case is discussed by Meerschaert et al. [57J . 
The coupled case is relevant as, in general, log-returns and waiting times are 
not independent [2Hj. Based on the results summarized in and discussed 
in 58 59 , it is possible to prove the following theorem for the coupled case: 



Theorem 



Let <p(£, t) be the (coupled) joint probability density of a CTRW. If, under the 
scaling £ — > h£ and r — > rr, the Fourier-Laplace transform of tp(£, r) behaves 
as follows: 

<?hA K > s ) = 0(hK,rs) (40) 
and if, for h — ► and r — * 0, the asymptotic relation holds: 

<Ph,ri K ' s ) = rs ) ~ 1 ~ mIM" - "(rs) 13 , (41) 

with < a < 2 and < < 1. Then, under the scaling relation /i/i" = vr$ , 
the solution of the (scaled) coupled CTRW master (integral) equation, eq. 0, 
Ph,r{x,t), weakly converges to the Green function of the fractional diffusion 
equation, u(x,t), for h — > and r — ► 0. 

Proo/ 

The Fourier-Laplace transform of the scaled conditional probability density 
p h , r (x,t) is given by: 

p fc , r («,a) = ~- -. (42) 

s 1 — ip[hn, rs) 

Replacing eq. (|41|l in eq. (|42|l and observing that ip( s ) = <?(0, s), one asymp- 
totically gets for small h and r: 

ur s 0-i 

which for vanishing h and r, under the hypotheses of the theorem, converges 
to: 

s P-i 

P0,0( K > S ) = = s p + ( 44 ) 

where u(k, s) is the Fourier-Laplace transform of the Green function of the 
fractional diffusion equation (see eq. (J2HJ)). The asymptotic equivalence in the 
space-time domain, between po t o(x, t) and u(x,t), the inverse Fourier-Laplace 
transform of u(k, s), is again ensured by the continuity theorem for sequences 
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of characteristic functions, after the application of the analogous theorem for 
sequences of Laplace transforms jSS]- There is convergence in law or weak 
convergence for the corresponding probability distributions and densities. 

An important consequence of the above theorem is the following corollary 
showing that in the case of marginal densities with finite first moment of wait- 
ing times and finite second moment of log-returns, the limiting density u{x, t) 
is the solution of the ordinary diffusion equation (and thus the limiting pro- 
cess is the Wiener process). The corollary can be used to justify the popular 
Geometric Brownian Motion model of stock prices, here with expected return 
set to zero. Again, in order to derive this result, no reference is necessary to 
the Efficient Market Hypothesis |551I55| . 

Corollary 

If the Fourier-Laplace transform of <p(£, r) is regular for k = and s = 0, and, 
moreover, the marginal waiting-time density, ipir), has finite first moment tq 
and the marginal jump density, A(£), is symmetric with finite second moment 
a 2 , then the limiting solution of the master (integral) equation for the coupled 
CTRW is the Green function of the ordinary diffusion equation. 

Proof 

Due to the hypothesis of regularity in the origin and to the properties of 
Fourier and Laplace transforms, we have that: 

\ / (o,o) V / (o,o) 

2 

= 1 - U — h 2 n 2 - T rs, (45) 

and as a consequence of the theorem, under the scaling h 2 a 2 /2 — Tor, one 
gets, for vanishing h and r: 

s + At 

corresponding to the Green function l|36|l for a = 2 and (3 = 1, that is the 
solution of the Cauchy problem for the ordinary diffusion equation. 

3 Summary and Outlook 

In this paper, a discussion of continuous-time random walks (CTRWs) has 
been presented as phenomenological models of tick-by-tick market data. 



Five Years of Continuous-time Random Walks in Econophysics 11 

Continuous-time random walks are rather general and they include compound 
Poisson processes as particular instances. Well-scaled limit theorems have been 
presented for a rather general class of CTRWs. 

It is the hope of this author that this paper will stimulate further research 
on high-frequency econometrics based on the concepts outlined above. There 
are several possible developments. 

First of all, one can abandon the hypothesis of i.i.d. log-returns and waiting 
times and consider various forms of dependence. In this case, it is no longer 
possible to exploit the nice properties of Laplace and Fourier transforms of 
convolutions, but, still, Monte Carlo simulations can provide hints on the 
behaviour of these processes in the diffusive limit. 

A second possible extension is to include volumes as a third stochastic 
variable. This extension is straightforward, starting from a three-valued joint 
probability density. 

A third desirable extension is to consider a multivariate rather than uni- 
variate model that includes correlations between time series. 

The present author is currently involved in these extensions and is eager 
to know progress in any direction by other independent research groups. He 
can be contacted at scalas@unipmn.it 
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